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Abstract 

We apply relativistic many-body methods to compute static differential polarizabilities for tran- 
sitions inside the ground-state hyperfine manifolds of monovalent atoms and ions. Knowing this 
transition polarizability is required in a number of high-precision experiments, such as microwave 
atomic clocks and searches for CP-violating permanent electric dipole moments. While the tra- 
ditional polarizability arises in the second-order of interaction with the externally-applied electric 
field, the differential polarizability involves additional contribution from the hyperfine interaction 
of atomic electrons with nuclear moments. We derive formulas for the scalar and tensor polarizabil- 
ities including contributions from magnetic dipole and electric quadrupole hyperfine interactions. 
Numerical results are presented for Al, Rb, Cs, Yb + , Hg + , and Fr. 

PACS numbers: 31.15.Ar,31.25.-v,32.60.+i 
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I. INTRODUCTION 



When an atom is placed in an external electric field, its energy levels shift due to the 
Stark effect. For states of definite parity, the effect arises in the second order in the inter- 
action of atomic electrons with the external E-field. The energy shift SE a is conventionally 
parameterized in terms of the polarizability of the atomic state a a , 

5E a = -\<si a £l , (1) 

where Sq is the strength of the applied E-field. 

The polarizability depends on atomic electric-dipole D matrix elements and energies E 

a ^M^f. (2) 



Eh — E a 

The sums are over a complete atomic eigen-set and the z-axis has been chosen along the 
E-field. On general grounds, we may decompose the polarizability for a state \nFM F ) of 
the total angular momentum F and its projection M F into the following contributions, 

s 3M F 2 - F(F + 1) T . . 

OinFMp - Ot nF H F(2F - 1) ^ ' ' ' 

Here the superscripts S and T distinguish the scalar and tensor parts of the polarizability. 
The "polarizabilities" a^ F and a^ F no longer depend on the magnetic quantum number M F . 

In this paper we focus on a difference of polarizabilities between two states nF' and nF 
of the same hyperfine manifold of states of total orbital angular momentum J = 1/2. Such 
calculations require additional care. Indeed, we are considering the Stark shift of hyperfine 
levels attached to the same electronic state. To the leading order, the shift is determined 
by the properties of the underlying electronic state. However, because the electronic state 
for both hyperfine levels is the same, the scalar Stark shift of both levels is the same. An 
apparent difference between the two levels is caused by the hyperfine interactionfHFI), and 
the rigorous analysis involves so-called HFI- mediated polarizabilities (see, e.g., [l]). Similar 
arguments hold for the tensor part of the polarizability. taken with its prefactor in 

Eq. ([2]), is an expectation value of an irreducible tensor operator of rank 2; it simply vanishes 
for J = 1/2 states due to the angular selection rules. Only the HFI coupling of nuclear and 
electronic momenta (F = I + J) leads to nonzero values of the tensor polarizability. 

Early works on Stark shifts of transition frequencies within hyperfine manifolds include 



Refs. 



2|-|6j. More recent interest to this problem was motivated by the Stark shifts of 



the hyperfine transition frequency due to the ambient black-body radiation (BBR) |7] . The 
BBR shift is a major systematic correction in microwave clocks, especially the 133 Cs primary 
frequency standard This motivated the most precise measurement of the DC Stark shift 
in a Cs fountain |9[. The relevant Stark shifts were a subject of many recent works (see, 



e.g., state-of-the-art calculations Ref. 10|, and references therein). 



Perhaps the most complete earlier theoretical treatment within the third-order (two 
electric-dipole couplings and one HFI) perturbation theory was given by Sandars 4j in 
1967. However, only recently (i.e., four decades later), a sign mistake in the expression of 
Ref. (4) for the tensor part of the HFI-mediated polarizability was discovered [ijl, Q] (a 
correct result for tensor polarizability of Tl was obtained earlier in Ref. |14j]). 

This sign error is directly relevant to extracting BBR correction from high-precision ex- 
periments. Notice that due to the isotropic nature of the BBR, the BBR clock shift is 
expressed in terms of the scalar part of the HFI-mediated polarizability. Moreover, charac- 
teristic frequencies of room-temperature BBR are well below excitation energies of atomic 
transitions thereby justifying replacing frequency-dependent polarizability with DC polar- 
izability [7j. Accordingly, the modern value of the BBR correction for the Cs clock is based 
on a measurement which was carried out in a DC E-field. However, the measured Stark 
shift involves a combination (j3]) of both scalar and tensor polarizabilities. Therefore to ar- 
rive at the BBR shift, one needs to remove the contribution due to the tensor polarizability. 
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13] becomes relevant. 



Clearly, the sign mistake discovered in 

Here we extend our earlier treatment of the HFI-mediated polarizabilities [l|, Q, with 
a specific focus on the tensor polarizabilities. Compared to Refs. 12, we employ a fully- 
relativistic formalism and evaluate tensor polarizabilities for several atoms and ions using 
modern relativistic many-body methods. We independently confirm that indeed, Ref. [4] 
had a sign mistake, requiring reinterpretation of measurements 9J. We also evaluate the 
tensor polarizability for the secondary frequency standard based on Rb atoms. 

Another motivation for our work comes from searches for the so far elusive permanent 
electric-dipole moments (EDM). Non- vanishing EDMs violate both time- and parity-reversal 
symmetries. Planned experiments will be carried out with Fr atoms [ija]. These atoms will 
be placed in a strong electric field and so-far unknown M^-dependent tensor polarizabilities 
would contribute to the error budget of the EDM search. Our computed values for Fr 
isotopes will aid the design and interpretation of these planned experiments. 



The paper is organized as follows. In Section [TT] we derive fully-relativistic third-order 
formulae for HFI-mediated tensor polarizabilities. Section IIHI presents details of numeri- 
cal evaluation within relativistic many-body theory. Finally, the results are discussed and 
compared with literature values in Section IIVI Unless specified otherwise, atomic units, 
m e = h = \e\ = 1 are used throughout. 



II. THEORETICAL SETUP 

We are interested in transitions between two hyperfine components of the same elec- 
tronic states. Below we employ the conventional labeling scheme for the atomic eigenstates, 
\n (U) FMp), where I is the nuclear spin, J is the electronic angular momentum, and F is 
the total angular momentum, F = J + 1. Mp is the projection of F on the quantization axis 
and n encompasses the remaining quantum numbers. 

As discussed in the introduction, computation of the transition polarizability for J = 1/2 
hyperfine-manifolds requires third-order analysis. This involves two perturbations due to 
the externally-applied electric field, Ve = — D ■ £ , and hyperfine interaction Uhfi- These 
perturbations may be chained into three distinct diagrams (see Fig. [TJ. Additionally, there 
is a residual (or normalization) diagram. 



— a 
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— < 
— < 
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(b) 



(c) 



(d) 



FIG. 1. Diagrams representing the complete third order hyperfine- mediated polarizability. Each 
diagram contains a hyperfine interaction (dotted line capped with a filled square) in addition to 
two interactions with the external electric field (dashed line capped with an empty triangle). The 
diagrams correspond to the (a) top, (b) center, (c) bottom, and (d) normalization terms discussed 
in the text. 



We would like to stress the importance of a consistent treatment of the HFI-mediated 
polarizabilities (i.e., including all the diagrams in Fig. [TJ. Consider a general expression for 
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the scalar polarizability, 

2 (nFMp | D k \i)(i \D k \nFM F ) 



k=x,y,z i=\mFiMi) 



Here all the involved states are the hyperfine states. While this requires that the energies 
include hyperfine splittings, it also means that the wave-functions incorporate HFI to all 
orders of perturbation theory. Including the experimentally-known hyperfine splittings in 
the summations is straightforward but limiting ourselves to this approximation would ex- 
clude the HFI corrections to the wave-functions. By expanding the energy denominators, 
we observe that including HFI into energies would only recover the residual diagram and 
partially the center diagram. We find that the remaining contributions are of the same 
order, and limiting computations to HFI-induced energy shifts only is hardly justified. 

Previously, we derived equations for dynamic HFI-mediated polarizabilities of hyperfine 
states in Ref. Clearly, static polarizabilities can be obtained by setting laser frequency 
to zero in the derived formulas. There is, however, one important addition to the formulae 
presented in Ref. [lj]: the HFI operator in that paper was truncated at the magnetic-dipole 
interaction. Here we additionally include the HFI coupling due to the electric-quadrupole 
nuclear moment. This contribution to tensor polarizabilities becomes increasingly important 
for heavier atoms. Note that the electric quadrupole contribution to scalar polarizabilities 
and thermal shift of states with total momentum j = 1/2 is zero. This can be explained 
in the following way. Scalar polarizability can be separated from total energy shift by 
averaging over directions of electric field. One cannot make non-zero scalar (energy shift is a 
scalar) from remaining electron angular momentum 1/2 and nuclear quadrupole, jajbQab = 
c a CbQab = (squared Pauli matrix sigma is reduced to delta symbol and antisymmetric 
linear tensor with a, Q a b is symmetric with zero trace, so that ^ a Q aa = 0). 

On general grounds, the rotationally-invariant hyperfine interaction between atomic elec- 
trons and nuclear moments may be written as a sum over scalar products of irreducible 



tensor operators (we follow notation of Ref. 16|) 



y HFI = ^Arw.r (7V) . (5) 



N 



Here the irreducible tensor operators J\f^ and act in the space of nuclear and elec- 
tronic coordinates respectively, with N being their ranks. The nuclear magnetic moment is 



conventionally denned as 

I2=(I,M I = I\^ 1) \I,M I = I) (6) 
and the nuclear electric-quadrupole moment as 

Q = 2(1, Mi = I\J\f^\I, M 7 = J) . (7) 

In the formulas below we require the reduced matrix element of the nuclear moment operator 
in the nuclear basis. For magnetic dipole, this is related to the nuclear magnetic g-factor as 

ai|A/- (1) ||/> = \y/{21) (2/ + l)(2/ + 2) gtm, 
H n being the nuclear magneton and \i = glfi n . For the electric-quadrupole moment, 

The components of relevant electronic tensors are 

i^2 (a • Cff (r 



T 



(i) 



where r is the electronic coordinate, a are the Dirac matrices, and (f ) and C\ (f ) are 
normalized vector spherical harmonic and normalized spherical harmonic functions, respec- 
tively hj|- 

n 

A derivation similar to Ref. [lj results in the scalar and the tensor polarizabilities given 
by (here [F] = 2F + 1) 

1 1 



(o) 

a nF ~ AT ri—T a nFi 



a nF 



(2F) (2F - 1) 1 1/2 



3(2F + l)(2F + 2) (2F + 3) 



(2) 



with the reduced polarizabilities a *Q being sums over values of individual diagrams of Fig. CEJ 



a 



(K) 
nF 



E( 2c Sv (T) + a%] N (C) + a% (O)) , (8) 



N 

where we used the equality of the top and bottom diagrams. 
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The angular reduction of individual diagrams leads to expressions 



«i^(T) = [F] V / M 



> J+Ja 



I I N 



JaJb \ J a J F 
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«S J v(C) = [F] A /M^^[J J ](-l) 2 ^ 
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• ,M ' i> 27 - J ) N 1 1 ^ ( n J\\T W \\nJ)(I\\Af^\\I)x 



<fn(0) = (-1Y 



F J J 
J J K 



) K J J l /V A ') 



1^11 



xO£'(nJ). 



F F I 

Finally, the universal (these are independent on F, i.e., the clock level) reduced sums are 
Tf)(nJN) = 2{I\\N {N) \\I) ^ nJ ^ riN) ! ( n a J a\\ D W n b J b) (n b J b \\D\\nJ) ^ 



n a ,n b =£n 



(E - E a ) (E - E b 



Cf) {nJ, N) = 2 || /) ^ (nJ| |£>| |w a J a ) (w a J a \ |T (jV) | K J h ) (n b J b \\D\\nJ) ^ 



n a ,n b ^n 



(E - E a ) (E - E b 



Ofj (n 



j^_ 2 Y^ ( n J\\ D \\ n aJa)( n aJa\\ D \\ n J) 

^ (E — E a ) 2 



(11) 



In these sums the values of the total orbital momenta of intermediate states J a and J b are 
fixed. 



III. NUMERICAL EVALUATION 



To perform the calculations we use an ah initio approach which has been described in 
detail in Ref . 18) . In this approach high accuracy is attained by including important many- 



body and relativistic effects. 
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Calculations start from the relativistic Hartree-Fock (RHF) method in the V N ~ X approx- 
imation. This means that the initial RHF procedure is done for a closed-shell atomic core 
with the valence electron removed. After that, the states of the external electron are calcu- 
lated in the field of the frozen core. Correlations are included by means of the correlation 



potential method [19|. We use the all-order correlation potential S for Rb, Cs, and Fr and 
second-order correlation potential E^ 2 ^ for Al, Yb + , and Hg + . The all-order £ includes two 
classes of the higher-order terms: screening of the Coulomb interaction and hole-particle 
interaction (see, e.g., (2^] for details). 

To calculate £ and X^ 2 -* we need a complete set of single-electron orbitals. We use the 

n 

B-spline technique [2JJ to construct the basis. The orbitals are built as linear combinations 
of 40 B-splines of order 9 in a cavity of radius 40ob. The coefficients are chosen from the 
condition that the orbitals are the eigenstates of the RHF Hamiltonian H Q of the closed- 
shell core. The all-order £ operator is calculated with the technique which combines solving 
equations for the Green functions (for the direct diagram) with the summation over complete 
set of states (exchange diagram) 20]. The second-order £( 2 ) operator is calculated using 
direct summation over complete set of states. 

The correlation potential £ is then used to build a new set of single-electron states, the 
so-called Brueckner orbitals. This set is to be used in the summation in equations ((9]), ( ITOl) 
and (TTTj) . Here again we use the B-spline technique to build the basis. The procedure is 
very similar to the construction of the RHF B-spline basis. The only difference is that new 
orbitals are now the eigenstates of the Ho + ^ Hamiltonian. 

Matrix elements of the HFI and electric dipole operators are found by means of the 



time-dependent Hartree-Fock (TDHF) method [19j, |22j . This method is equivalent to the 
well-known random-phase approximation (RPA). In the TDHF method, the single-electron 
wave functions are presented in the form ip = ip + Sip, where ipo is the unperturbed wave 
function. It is an eigenstate of the RHF Hamiltonian Hq\ (Ho ~ ecOV'o — 0. 5ip is the 
correction due to external field. It can be found be solving the TDHF equation 

(H - e )5ip = -5e% - i^o - 5V N ~^ , (12) 

where 5e is the correction to the energy due to external field (5e = for the electric dipole 
operator), F is the operator of the external field (Vhfi or — D£), and ^V^ -1 is the correction 
to the self-consistent potential of the core due to external field. 



S 



The TDHF equations are solved self-consistently for all states in the core. Then the 
matrix elements between any (core or valence) states n and m are given by 

(iinlF + SV*- 1 ^). (13) 

The best results are achieved when ip n and ip m are the Brueckner orbitals computed with 
the correlation potential E. 

We use equation ( II 3p for all HFI and electric dipole matrix elements in evaluating the top, 
bottom, and center diagrams (Eqs. (!9|). (|T0|) .( fTT]) ) except for the ground state HFI matrix 
element in the normalization diagram where we use experimental data. The results are 
presented in section [TV] 

IV. RESULTS AND DISCUSSION 

Table [J shows the results of calculation of scalar and tensor hyperfme mediated polariz- 
abilities for two hyperfme components of the ground state of some isotopes of Al, Rb, Cs, 
Yb + , Hg + , and Fr. Magnetic dipole and electric quadrupole HFI contributions to the tensor 
polarizabilities (a T ) are shown separately. Notice that there is no electric quadrupole contri- 
bution to the scalar polarizabilities (a 5 ), as discussed previously. Another interesting thing 
to note is that electric quadrupole does not contribute to the frequency shift of the hyperfme 
transition. Although the tensor polarizabilities are different for states with F = I + 1/2 and 
F — I — 1/2, the energy shifts, which includes F-dependent factors (see Eq. Q) are the 
same. This is only true for states with Mp = 0, where Mp is projection of F. 

The accuracy of the calculations is different for scalar and tensor polarizabilities, being 
a few per cent for scalar polarizabilities and about 30% for tensor polarizabilities. This is 
because we include only Brueckner-type correlations, or correlations which can be reduced 
to redefinition of single-electron orbitals. Such approximation works very well for s\/2 and 
Pi/2 states which dominate in the scalar polarizabilities. For tensor polarizabilities large 
contribution comes from p 3 / 2 and d 3 / 2 states where Brueckner approximation is not so good, 
especially for the HFI matrix elements. To get accurate results for tensor polarizabilities 
one has to include structure radiation and other non-Brueckner higher-order correlation 
corrections. This goes beyond the scope of present work. 

To illustrate the accuracy of the calculations we would like to compare our results to 
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TABLE I. Third-order hyperfine static polarizabilities of single-valence atoms. 



z 


Atom 


.4 


I 




Ql 

[b] 


F 


a s 


io- 


^ 10 Hz/(V/m) 2 


a Te _ 


13 


Al 


27 


5/2 


3.6415 


0.14 


3 


0.0158 


0.2683 


-0.0119 


0.2563 
















2 


-0.0222 


-0.1073 


-0.0096 


-0.1169 




37 


Rb 


85 


5/2 


1.3530 


0.27 


3 


0.4599 


-0.0073 


-0.0070 


-0.0143 
















2 


-0.6439 


0.0029 


-0.0056 


-0.0027 




37 


Rb 


87 


3/2 


2.7510 


0.132 


2 


0.9332 


-0.0148 


-0.0034 


-0.0182 


-0.0234 














1 


-1.5554 


0.0025 


-0.0017 


0.0007 


0.0009 


55 


Cs 


133 


7/2 


2.5820 


-0.004 


4 


1.9770 


-0.0262 


0.0002 


-0.0260 


-0.034 














3 


-2.5419 


0.0141 


0.0002 


0.0143 


0.0184 


70 


Yb+ 


171 


1/2 


0.4940 


1 


1 


0.0844 


-0.0023 





-0.0023 



















-0.2533 













70 


Yb+ 


173 


5/2 


-0.6775 


2.8 


3 


-0.1158 


0.0031 


-0.0390 


-0.0359 
















2 


0.1621 


-0.0012 


-0.0312 


-0.0325 




80 


Hg+ 


199 


1/2 


-0.5603 


0.4 


1 


0.0271 


-0.0018 





-0.0018 



















-0.0814 













80 


Hg+ 


201 


3/2 


0.5060 


0.8 


2 


-0.0300 


0.0020 


-0.0014 


0.0005 
















1 


0.0501 


-0.0003 


-0.0007 


-0.0011 




87 


Fr 


211 


9/2 


4.0005 


-0.19 


5 


6.9771 


-0.1459 


0.0188 


-0.1271 


-0.1633 














4 


-8.5275 


0.0908 


0.0175 


0.1083 


0.1392 


87 


Fr 


221 


5/2 


1.5800 


-0.98 


3 


2.7556 


-0.0576 


0.0970 


0.0394 


0.0506 














2 


-3.8579 


0.0230 


0.0776 


0.1006 


0.129 


87 


Fr 


223 


3/2 


1.1700 


1.17 


2 


2.0406 


-0.0427 


-0.1158 


-0.158 


-0.204 














1 


-3.4010 


0.0071 


-0.0579 


-0.0508 


-0.0653 



a Reference 
b Magnetic dipole HFI contribution 
c Electric quadrupole HFI contribution 
d Total tensor polarizability, a T — + 
c Corrected tensor polarizability (see text for discussion) 
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TABLE II. Stark frequency shift: comparison with previous calculations. 



Atom 


This work 


k, [10- 10 Hz/(V/m) 2 l 

' L 1 \ 1 / 1 


Ref. [18] 


87 Rb 


-1.24(1) 




-1.24(1) 


133 Cs 


-2.26(2) 




-2.26(2) 


171 Yb 


-0.167(9) 




-0.171(9) 


199 Hg 


-0.056(3) 




-0.060(3) 



previous calculations and available experimental data. There were numerous calculations 
and measurements of the Stark frequency shift for the hyperfine transitions of the ground 



state of Cs, Rb, and other atoms and ions (see, e.g., Ref. 18[ and references therein). The 



results are usually presented in terms of the coefficient k related to frequency shift by 

5AE = k£l , (14) 

where AE is the energy interval between two hyperfine components of the ground state 
and 5AE is its change due to electric field Sq. One can find the values of k for all atoms 
and ions from Table [I] using the polarizabilities from this table and expressions (PO) and (El). 
Corresponding values for 87 Rb, 133 Cs, 171 Yb + , and 199 Hg + are presented in Table [III and 
compared to previous calculations of Ref. jl8| . Since both calculations are performed with 
the same method, we assume the same uncertainty. Some differences in central values for 
Yb and Hg are due to differences in the details of the calculations. These differences are 
within the declared uncertainty. 

As has been discussed above the accuracy for the tensor polarizabilities is lower. Our 
value of the tensor polarizability for the F = 4 state of cesium is —2.60 Hz(V/cm) 2 . This 



is about 30% smaller that the experimental value —3.34(2) (25) Hz(V/cm) 2 24j and the 



semiempirical value —3.72(25) Hz(V/cm) 2 13j. The most likely reason for the difference 



is the contribution from the non-Brueckner correlations which are not included in present 
work. Therefore, we assume 30% uncertainty for the calculated tensor polarizabilities in 
present work. 

To improve the predicted values for the tensor polarizabilities of Rb and Fr, which both 
have electron structure similar to those of cesium, we multiply the ab initio results for these 
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atoms by the factor 1.28, which is the ratio of the experimental tensor polarizability for Cs 
to the theoretical one. The resulting values are presented in the last column of Table HI 
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